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ABSTRACT 

We use new and archival Chandra data to investigate the X-ray emission from a 
large sample of compact hotspots of FRII radio galaxies and quasars from the 3C cat- 
alogue. We find that only the most luminous hotspots tend to be in good agreement 
with the predictions of a synchrotron self-Compton model with equipartition mag- 
netic fields. At low hotspot luminosities inverse-Compton predictions are routinely 
exceeded by several orders of magnitude, but this is never seen in more luminous 
hotspots. We argue that an additional synchrotron component of the X-ray emission is 
present in low-luminosity hotspots, and that the hotspot luminosity controls the ability 
of a given hotspot to produce synchrotron X-rays, probably by determining the high- 
energy cutoff of the electron energy spectrum. It remains plausible that all hotspots 
are close to the equipartition condition. 

Subject headings: galaxies: active - X-rays: galaxies - X-rays: quasars - radiation 
mechanisms: non-thermal 



1. Introduction 

The hotspots of powerful FRII radio sources, as observed in the radio, have long been believed 
to be the observable consequence of a strong terminal shock at the end of the relativistic jets that 
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feed the radio lobes. This picture has to be modified somewhat in the Ught of the fact that a large 
number of lobes contain more than one hotspot when observed at high resolution (Laing 1982; 
Leahy et al. 1997; Hardcastle et al. 1997). Conventionally the most compact feature, which is 
universally the one that lies at the end of the jet when the jet termination can be observed, is 
called the 'primary' hotspot, while the others are known as 'secondary' hotspots. Since standard 
shock acceleration and energy loss models are often good fits to the radio-to-optical spectra of both 
primary and secondary hotspots (e.g., Meisenheimer et al. 1989), it seems likely that in many cases 
the secondary hotspots are also being powered by bulk kinetic energy from the jet. 

For some time it has been clear that it is not possible to explain the X-ray emission from 
hotspots with a single model. Some objects show hotspot X-ray emission with a spectrum consis- 
tent with the predictions of a synchrotron self-Compton (SSC) model, in which the synchrotron- 
emitting electrons inverse-Compton (IC) scatter synchrotron photons into the X-ray band; in these 
objects there is good agreement between the observed flux density and the predictions of an SSC 
model with a magnetic field close to the value expected for equipartition of energy between the 
magnetic field and the radiating electrons, and so their X-ray emission is attributed to the SSC pro- 
cess (e.g. Harris, Carilli & Perley 1994; Hardcastle et al. 2002, hereafter H02). On the other hand, 
objects such as 3C 390.3 (Harris, Leighly & Leahy 1998) and Pictor A (Roser & Meisenheimer 
1987) show X-ray emission that is clearly much stronger than the SSC model would predict if the 
magnetic field had the equipartition value, together with a spectrum that is steeper than the low- 
frequency radio spectrum, and hence too steep to be inverse-Compton (Wilson, Young & Shopbell 
2001; H02). In some, but not all, of these cases, a simple synchrotron spectrum (by which we 
mean a single power law or a broken, steepening power law in frequency) is a good fit to the ra- 
dio, optical and X-ray data points. In addition, there are several sources (the best example being 
3C 351: H02) where the X-ray structure is clearly different from that seen in the radio maps, which 
is impossible in a simple SSC model with a homogeneous magnetic field and electron distribution. 
In a synchrotron model for some or all of the X-rays, differences in the spatial structure are to be 
expected, since the synchrotron loss timescale for X-ray emitting electrons (tens of years in a typ- 
ical equipartition magnetic field) is orders of magnitude less than that for radio-emitting electrons 
(> 10^ years); in fact, in a non-steady-state situation (as expected from numerical simulations, e.g. 
Tregillis, Jones & Ryu 2001) time-varying differences in both spatial and spectral distributions of 
the radio and X-ray emitting electrons are more or less required by the physics. 

Until now, however, it has not been clear why some hotspots' X-ray emission is adequately 
modeled by the SSC process with an equipartition field, while others require an additional com- 
ponent or a lower than equipartion field strength. Suggested explanations have involved (1) a 
lower magnetic field strength in the X-ray bright hotspots, which both increases IC emission and 
increases the loss lifetime of X-ray synchrotron-emitting electrons (Brunetti et al. 2001); or (2) 
the effects of differential relativistic beaming, due to decelerating bulk motions in the hotspots. 
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on the synchrotron and IC spectra, in particular the fact that fast-moving parts of the flow see the 
slow-moving downstream flow Doppler-boosted (Georganopoulos & Kazanas 2003), which helps 
to account for the fact that many of the early detections of X-ray bright hotspots involved broad- 
line radio galaxies or quasars (H02), which in unified models should lie relatively close to the line 
of sight. The small number of published detections of X-ray hotspots has made it difficult to arrive 
at a definitive answer. 

Since the appearance of earlier work aimed explicitly at detecting hotspots (Harris et al. 2000; 
Hardcastle et al. 2001a, 2002; Brunetti et al. 2002) a number of hotspots have been detected in 
Chandra observations of FRII radio galaxies and quasars made for other purposes (e.g., Donahue 
et al. 2003; Crawford & Fabian 2003). This has motivated us to collate all existing data on the 
hotspots of 3C radio sources from the Chandra archive and analyse them in a systematic way, with 
the aim of determining trends and testing models. In this paper we report our results. 

Throughout the paper we use a cosmology in which Hq = 70 km s^ Mpc~^, = 0.3 and 
Qa = 0.7. Spectral indices a are the energy indices and are defined in the sense S^, oc 

2. Data and analysis 

We searched the public Chandra archives for all observations of FRII radio sources in the 
3C catalogue made with the CCD Imaging Spectrometer (ACIS), supplementing them with a few 
observations that we have access to and are not yet public. To make our sample as large as possible, 
we included any 3C FRII source listed in the archives; this means that we made use of several 
sources that are not in the better-defined 3CR (Spinrad et al. 1985) or 3CRR (Laing et al. 1983) 
samples. We restricted ourselves to the 3C sources simply because they almost all have good 
radio observations available in the public NRAO Very Large Array (VLA) archives. We excluded 
compact steep-spectrum sources whose hotspots would not be resolved from the AGN or lobe 
emission with Chandra. These selection criteria gave us a sample of 36 sources (Table 1). We next 
obtained electronic radio maps for all the sources in our sample. Where we did not already have 
access to a good, high-resolution radio map, or to appropriate published radio flux densities, we 
retrieved the best available data from the VLA archive. In selecting the observations we preferred 
data at 5 and 8 GHz in the A configuration (the largest configuration of the VLA), which gives 
good sensitivity, good separation of lobes and hotspots due to their spectral index differences, and 
an angular resolution comparable to that of Chandra. Properties of the radio data used and maps 
made are given in Table 2. 

We then used the radio maps as a guide to search for emission from hotspots. We defined 
a hotspot less strictly than some other workers (e.g.. Bridle et al. 1994) both in order to make 
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the comparison between radio and X-ray simpler and because our aim is to include the structure 
apparently associated with the jet termination whenever possible. In practice, we considered any 
relatively compact, isolated radio feature that was significantly brighter than its surroundings to 
be a hotspot (cf. Leahy et al. 1997). However, we excluded any emission that we considered to 
be associated with a jet; any compact X-ray feature positionally coincident with, or closer to the 
nucleus than, a linear radio feature that met the definition of Bridle & Perley (1984), was considered 
to be a 'jet knot' rather than a hotspot. By doing this we hoped to select only features associated 
with the termination of the jet, and to avoid effects thought to be due to highly relativistic bulk 
motions, as seen in the X-ray jets of some quasars (e.g., Tavecchio et al. 2000). Later in the paper 
(§5) we shall return to the question of whether hotspots tell us anything about the emission from 
jet knots or jets in general. 

Using the radio data as a guide, we were able to identify a number of previously unreported 
X-ray counterparts to hotspots, in the 10 FRIIs 3C6.1, 3C47, 3C 109, 3C 173.1, 3C228, 3C 321, 
3C324, 3C334, 3C403 and 3C452. The resuhs on 3C403 and 3C228, 334 will be reported in 
more detail elsewhere (respectively in R. Kraft et al., and D.M. Worrall et al., in preparation); 
images of the newly detected X-ray hotspots for the other sources are presented in Appendix A. 
In almost all cases there were a few tens of total counts in the hotspots in the 0.5-5 keV range. 
Since this is too few to fit spectra, we used the standard Chandra analysis software CIAO to gen- 
erate redistribution matrix (RMF) and ancilliary response (ARF) files appropriate for the hotspots, 
using the PSEXTRACT tool, corrected for the time-dependent excess ACIS absorption using ap- 
ply _acisabs, and then used the model-fitting software XSPEC to determine the normalization of 
a power law, with a — 0.5 and Galactic absorption, that reproduced the observed net count rate. 
We chose this power-law index because our aim was to test the validity of the SSC model, which 
'predicts' a = 0.5, on the basis of the assumption that the low-energy electron energy index has 
the value 2.0 associated with particle acceleration at a non-relativistic strong shock (as appears to 
be the case in some, though not all, well-studied hotspots: Meisenheimer, Yates & Roser 1997). 
However, the choice of a makes relatively little difference to the normalization of the power law, 
and thus to the inferred 1-keV flux density in the observer's frame; using a — 1 would increase 
the inferred flux density by between 10 and 20%. The 1-keV flux densities assuming a — 0.5 for 
each source are given in Table 3. 

For the 10 sources 3C9, 3C 184, 3C200, 3C212, 3C215, 3C219, 3C 220.1, 3C401, 3C 427.1 
and 3C 438 we found no X-ray emission associated with the hotspots, as previously reported in 
some cases. Where there was a compact (arcsec or sub-arcsec) radio hotspot present we determined 
3(7 upper limits on the corresponding X-ray flux density in a 10-pixel detection cell based on 
Poisson statistics and the local background count rate. Some of these sources (e.g. 3C401 and 
3C438, Hardcastle et al. 1997) exhibit no compact hotspots in the radio, and we elected not to 
determine an upper limit on their emission, as the selection of an appropriate X-ray region is 
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difficult. For the same reason, we did not determine an upper limit for the barely resolved source 
3C 1 84, which shows X-ray emission coincident with the radio lobe but not particularly the hotspot 
(Belsole et al. 2004). 

We re-examined previously reported hotspot detections in the 16 sources 3C 123 (Hardcastle 
et al. 2001a), 3C 179 (Sambruna et al. 2002), 3C207 (Brunetti et al. 2002), 3C254 (Donahue et 
al. 2003), 3C 263 (H02), 3C 265 (Bondi et al. 2003), 3C 275.1 (Crawford & Fabian 2003), 3C 280 
(Donahue et al. 2003), 3C281 (Crawford & Fabian 2003), 3C294 (Fabian et al. 2003), 3C295 
(Harris et al. 2000), 3C303 (Kataoka et al. 2003), 3C330 (H02), 3C351 (Brunetti et al. 2001; 
H02), 3C 390.3 (Harris et al. 1998), and Cygnus A (3C405: Harris et al. 1994). In all but one 
case we confirmed the existence of one or more compact X-ray features associated with the radio 
hotspots (the exception is 3C281, where the previously reported X-ray emission appears to be 
diffuse and associated with the lobe, and is most likely due to IC scattering of cosmic microwave 
background (CMB) photons by the lobe rather than SSC from the hotspots). Where a 1-keV flux 
density had been previously determined from spectral fitting, we make use of that in Table 3. Oth- 
erwise, we adopted the same procedure as for the newly detected sources described above. Finally, 
for all the newly detected and known sources, we determined upper limits, again as described 
above, for any compact hotspots that were not detected (e.g. in the lobe on the opposite side of 
the nucleus to the known hotspot). In the process of doing this we found one additional hotspot, 
3C 123 W, that was formally significantly detected; although the situation in this source is confused 
by the presence of strong, unrelaxed cluster emission (Hardcastle et al. 2001a) we added it to the 
sample as a detection, for consistency with the other sources. Because of its intrinsic interest and 
the extreme nature of its X-ray hotspot, we added Pictor A (Wilson, Young & Shopbell 2001) to 
the sample (it is not in 3C because of its low declination, but meets the other selection criteria). 
The overall final sample thus contains 37 sources (Table 1) with 65 X-ray hotspot flux densities or 
upper limits. All the fluxes and upper limits are tabulated in Table 3. 

Finally, for the sources with detected X-ray hotspots, we used the Hubble Space Telescope 
(HST) archive to search for optical counterparts. Optical emission is important because it con- 
strains the spectrum between radio and X-ray; early work on SSC hotspots was supported by the 
observation that a one-zone synchrotron model could not be fitted through the radio and X-ray data 
points because of the optical constraints (e.g., Harris et al. 1994). We identified two new candidate 
optical hotspot counterparts (in 3C 228 and 3C 275.1) and measured flux densities or upper limits 
for a number of other sources, using the IRAF package synphot to calculate the conversion factor 
between observed counts and flux density. Optical flux densities and frequencies are tabulated in 
Table 4. Sources where there were no archival HST observations, where the hotspot did not lie 
on the WFPC-2 CCDs, or where observational constraints such as a nearby bright star or cosmic 
ray contamination prevented us from obtaining a flux density, are not tabulated. We also tabulate 
a number of flux densities, largely based on ground-based observations, taken from other papers. 
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either in the literature or in preparation. 



3. Modeling and results 

The large number of detected X-ray hotspots is interesting in itself, given that calculations 
based on SSC emission at equipartition suggested that only the few brightest hotspots would be 
detected with Chandra (e.g., Hardcastle 2001). In order to assess quantitatively the extent to which 
the new detections conflict with an SSC model, we decided to fit a simple SSC model to all the 
hotspots and determine the ratio between the observed and predicted flux densities. We carried 
out this calculation using the code of Hardcastle, Birkinshaw & Worrall (1998); a brief sketch of 
the operation of this code is given in Appendix B. The code assumes a spherically symmetric, 
homogenous hotspot with an electron energy spectrum that can be described as a power law or 
broken power law. To determine the radius of the hotspot, we therefore fitted models consisting 
of a homogeneous sphere convolved with the restoring beam to the highest-resolution radio data 
available, in the manner described by H02. Where multi-frequency radio data were available for the 
hotspot, which was only true in the best-studied cases, we used them to fit a two-component power- 
law model with an energy spectral break corresponding to Act = 0.5 (Heavens & Meisenheimer 
1987); otherwise we assumed a single power law with a = 0.5 extending from the radio into the 
mm-wave regime. Unless good low-frequency radio constraints were available, we assumed that 
the minimum Lorentz factor of the electrons, 'jmin, was 1000; 7max was chosen to ensure that there 
was no spectral cutoff before the mm-wave region. These choices reflect what has been found 
in the best-studied hotspots, but are clearly no substitute for good, multi-frequency observations, 
particularly at high frequencies. However, we estimate that these choices make a difference at the 
level of at most ~ 10-20% (except in the rare cases where a ~ 10 GHz spectral cutoff is present). 
The results are particularly insensitive to the choice of 7inin> since decreasing this has two effects 
which act in opposite directions: more high-energy photons are scattered by the large additional 
population of low-energy electrons, but the overall electron energy power law normalization is 
reduced to maintain equipartition. Detailed spatial modeling, where high-resolution observations 
make it possible, also changes the results of SSC calculations at the 10-20% level (H02), so that 
overall the calculated value should be a good estimate of the true inverse-Compton prediction. 
The equipartition flux density prediction (taking into account both SSC and IC scattering of CMB 
photons) and the ratio R between the observed and predicted flux densities is tabulated in Table 3 
for each hotspot. Note that in almost all cases the flux due to the SSC process dominates over that 
due to IC scattering of the CMB (assuming no relativistic beaming) by an order of magnitude or 
more. For simplicity we shall often refer to the calculated fluxes as SSC fluxes in what follows. 

The tabulated values of R are calculated assuming that there are no protons, so that equipar- 
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tition is between the radiating electrons and magnetic fields only, and that the filling factor of the 
hotspots is unity. If we were to include an energetically dominant population of protons in equipar- 
tition in our model, it would reduce the predicted inverse-Compton emission (since the number of 
electrons decreases) and so increase the ratio R, possibly by a large factor. If the ratio of proton 
to electron energy densities were in the ratio of their rest masses, R values would increase by 
about a factor 30-70. Even if the proton to electron energy densities were of the order of their 
number ratios as observed in cosmic rays at the Earth (an energy-dependent factor of ~ 50-100, 
Longair 1992) we would expect R to increase by a factor ~ 5-10. A population of protons with 
the same total energy as the electrons has a less dramatic effect, increasing R by only about 40%. 
In any case, it is clear that introducing protons cannot solve the problem of high-i? hotspots. The 
predicted SSC inverse-Compton emission can be increased, and R can be decreased, if we have 
overestimated the volume or the filling factor, though (depending on the space-filling fluid) the ac- 
tual results of a low filling factor can be very geometry-dependent - if the electrons are confined to 
thin sheets, so that the probability of scattering is comparatively low, then the effects of low filling 
factor can be less than expected. Roughly (see Appendix B), to reduce i? by a factor of 1000, and 
so to make the most extreme observed hotspots consistent with being inverse-Compton emission 
at equipartition, we would need to reduce the volume or the filling factor by a factor ~ 10^^, and 
this neglects geometrical effects. Such low filling factors are clearly implausible. 

Fig. 1 shows the distribution of the ratio R and upper limits on R for the sample. Two points 
are immediately obvious: firstly, the detected X-ray flux density lies significantly above the IC 
prediction in most sources; secondly, there are few detected sources with R < 1, and there is a 
clear change in the distribution of sources at around R — 1. If the upper limits on X-ray flux for 
the non-detected sources all lie a long way above the true values, then this could change, but if it 
does not, the special status of i? = 1 implies that few sources have less X-ray emission than would 
be expected on the equipartition SSC/IC model. This could suggest either that all sources have 
SSC/IC X-ray emission at a level consistent with the equipartition prediction, together with some 
additional source of X-ray emission, or that there are departures from equipartition — some quite 
large — but that these are always in the sense that B < i?eq. The fact that we do not see many 
sources with R <^ 1 suggests that there are few or no hotspots with B ^ Beq, though without 
detections of all the hotspots we cannot be more definite. 

If the model we have used to predict the level of inverse-Compton emission is incorrect, then 
the special status of R — 1 would have to be a coincidence. This, as we have argued before 
(e.g. H02) gives us a reason to disfavor models with an energetically dominant proton population, 
or with consistently very low filling factor. Thus, for example, if the proton-to-electron energy 
density ratio were ~ 100, we would expect a source in equipartition, and emitting in X-rays only 
via the SSC process, to have R ^ 0.1-0.2 given our model (since we would be overpredicting the 
inverse-Compton emission); we would not see a special status for = 1 unless other parameters 
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(such as filling factor) conspired systematically to increase the value of R, which is inherently 
improbable. We emphasise that, given the small numbers and the width of the distribution of R 
values, this does not rule out moderate filling factors or a proton population within a factor of ~ 10 
of the energy density in the other components. 

What emission mechanisms are possible for the detected hotspots? We began by calculating 
another parameter, R', the ratio between the observed X-ray flux density and the flux predicted 
from a simple power-law extrapolation (with a^x = 10) from the radio data. Hotspots with both 
R^ I and R' ^ 1 would represent a problem for both synchrotron and inverse-Compton models. 
However, we find that R' is almost always less than 1, so there are no sources whose X-ray flux 
is impossible to explain with a synchrotron model in this sense. The plot of R' against R (Fig. 2) 
shows that even the most extreme X-ray hotspots, in terms of R, can readily be accounted for with a 
synchrotron model. There is a smooth distribution in parameter space, with no obvious bimodality, 
between sources with R' — 1, R ^ 1 (where a synchrotron model is natural) and R — 1, R' <^ 1 
(where an inverse-Compton model has tended to be adopted in earlier work). Insets in Fig. 2 show 
the very different SEDs of sources at the extreme ends of the distribution, and also illustrate the 
importance of optical constraints in determining the X-ray emission mechanism. To investigate this 
further, we used the available optical data or upper limits for X-ray detected hotspots to constrain 
their spectral shape. We calculated the quantities and aox, the two-point radio to optical and 
optical to X-ray spectral indices, for all the X-ray detected hotspots with optical flux densities or 
upper limits. A hotspot in which the optical to X-ray spectrum is flatter than the radio to optical 
spectrum (aox < ckro) cannot be described by a simple one-zone synchrotron model in which the 
spectrum steepens with increasing frequency. The difference between the two spectral indices is 
plotted in Fig. 3 as a function of R. It can be seen that the more extreme hotspots (large R values) 
all have «ox — ono > 0, and so are consistent with a synchrotron model. A non-synchrotron model 
is required by the optical data only for a few low-i? objects, where inverse-Compton emission is the 
accepted and most plausible X-ray mechanism. The existing data do not rule out a model in which 
synchrotron X-ray emission is important in a significant number of our target objects, although this 
is far from conclusive given the large number of optical non-detections. 

What determines the value of R for a particular hotspot? We noted that the early detections 
of SSC emission, such as Cygnus A, 3C295 and 3C 123, were all in luminous sources, while 
well-studied problematic sources such as 3C 390.3 and Pictor A are much lower in overall radio 
luminosity. Accordingly, we looked for a relationship between it! and total 178-MHz luminosity 
from the original 3C measurements or the revised values of Laing et al. (1983) (correcting by a 
factor 1.09 so as to bring the flux densities on to the scale of Baars et al. (1977), and using low- 
frequency spectral indices to correct to the rest frame), obtaining the plot shown in Fig. 4. The 
inverse correlation seen here appears to indicate a role of the source luminosity in determining R. 
The correlation is improved if we plot the luminosity of only the hotspot against R (Fig. 5); here 
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we have used measurements from the radio maps, correcting to a rest-frame frequency of 5 GHz 
by assuming a radio spectral index of 0.5. The improvement suggests that the relationship with 
hotspot luminosity is primary, and that the correlation with overall source luminosity arises because 
of the correlation between hotspot and source luminosity. It is important to realise that these plots 
are not necessarily an indication of a one-to-one correlation between R and hotspot or total source 
luminosity. Firstly, Chandra's sensitivity (around 0.1 nJy at 1 keV for the exposure times used in 
these observations) means that we would not expect to detect the SSC emission from the lowest- 
luminosity hotspots, so that it is observationally impossible to populate the bottom left-hand comer 
of Fig. 5, as shown by the dotted lines illustrating the observational limits. Secondly, there is a 
positive correlation between the predicted SSC flux density and the hotspot luminosity, since the 
sample is flux-limited and the SSC luminosity is a non-linear function of the hotspot luminosity (for 
a given hotspot size) and this increases the strength of the apparent correlation. However, there is at 
least one key result from this analysis; there are no hotspots with high luminosity and high R — we 
would certainly have been able to detect such hotspots if they existed. By contrast, low-luminosity 
hotspots appear to be able to have extremely high R values, though we cannot say definitely that 
all of them do. We find no other relationships between R and hotspot or source parameters such 
as hotspot angular or linear size, source size, redshift, or radio spectral index. However, there are 
relationships between R and other derived quantities such as equipartition magnetic field energy 
density and photon energy density (Fig. 6); these are not surprising, since all of the quantities are 
related to radio luminosity. We return to the possible physical significance of these relationships 
below (§4.2). 

Finally, we investigated the role of beaming in determining the X-ray brightness of hotspots by 
plotting the R parameter against the core prominence, defined here as the ratio of 5-GHz core flux 
density to (rest-frame) 178-MHz total source flux density. Core prominence is often used as a proxy 
of beaming (e.g. Orr & Browne 1982; Kapahi & Murphy 1990; Morganti et al. 1997; Hardcastle 
et al. 1999), relying on the assumption that the intrinsic fraction of the radio source flux emitted 
by the core is similar in all sources and that the observed variation in core prominence arises from 
relativistic beaming in the parsec-scale jet. Using the core flux densities tabulated in Table 1, we 
produced the plot shown in Fig. 7. This figure certainly shows a trend, in the sense that (as was 
already clear) many of the X-ray over-bright hotspots are in beamed sources, and often on the same 
side as a known one-sided radio jet. At the same time, there is clearly a good deal of scatter in any 
correlation — up to 2.5 orders of magnitude separate sources with similar core prominences — 
and there are sources that do not fit it at all, such as the low core-prominence, narrow-line source 
3C403 {R > 1000). If we plot core prominence against hotspot radio luminosity (Fig. 8) we see 
that there is a tendency for sources with high core prominences to have low-luminosity hotspots, a 
trend that can be explained entirely in terms of a bias towards broad-line radio galaxies and quasars 
at low redshifts in the parent sample, so that it is not clear that the trend seen in Fig. 7 is meaningful. 



-10- 



Although the quasars and broad-line radio galaxies (in unified models, the sources that should be 

most strongly affected by beaming) tend to lie at the upper edge of the envelope of R values for 
a given luminosity range, the dominant effect is the luminosity dependence. A partial Kendall's r 
analysis taking into account the upper limits, performed according to the prescription of Akritas & 
Siebert (1996) and using their code, shows that the correlation between R and core prominence is 
not significant at the 95% confidence level if the luminosity correlation is taken into account, while 
the correlation between R and luminosity is significant even given the core prominence relation. 
So, with the current data, we have no significant evidence for a relationship between hotspot R 
value and beaming, and the apparent correlation of Fig. 7 must be regarded as suggestive at best. 

4. Discussion 

What models can account for these observations? We examine several in turn. 



4.1. SSC with luminosity-dependent departures from equipartition 

One obvious possibility is that the high R values reflect a significant departure from equipar- 
tition; the belief that equipartition fields exist in hotspots is, after all, a result of the study of the 
most luminous hotspots (chosen because of their high predicted SSC flux densities). The departure 
from equipartition in terms of the ratio of equipartition to true magnetic field strengths, Bcq/ B, is 
approximately , so that the magnetic field strength would have to be a factor ~ 100 lower than 
the equipartition value in the hotspots with the highest R values. 

We can rule out a simple and attractive model in which such low magnetic fields account 
entirely for the luminosity-/? correlation. In this toy model, all hotspots have similar numbers of 
electrons and sizes. The radio luminosity Lr goes as 5^+", and the IC luminosity scales in the 
same way (so long as the synchrotron photon field remains dominant). But the equipartition pre- 
diction for SSC decreases more rapidly, since this depends on the equipartition estimate of electron 
density, which goes approximately as L^^^ (Appendix B), as well as linearly on the observed pho- 
ton density. So we should find that R oc L~^/^, which is not far from the observed slope. We should 
also find that the hotspot radio luminosity scales approximately as B^^, where B^eq is the magnetic 
field strength required to produce the observed X-ray emission by inverse-Compton processes, and 
this is also just about consistent with the data for our sources. However, we would not expect to see 
the observed correlation between 178-MHz total flux density and R (Fig. 4) in this picture, unless 
the same ratio between the true and equipartition magnetic fields persisted throughout the source. 
Very low fields in lobes are incompatible with observations of lobe inverse-Compton emission. 
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among other things; we see high-i? hotspots in sources whose lobes clearly do not show the same 
ratio between the observed and predicted inverse-Compton emission from CMB photon scattering 
(e.g. 3C403, R. Kraft et al. in preparation). 

More generally, there are several arguments that seem to us to disfavor an SSC model (to be 
accurate, SSC plus inverse-Compton scattering of the CMB) with a larger departure from equipar- 
tition in lower-power sources: 



• We know (see §1) that some high-i? hotspots' X-ray spectra and/or spatial properties are 
inconsistent with a pure SSC model (the spectra of the best-studied high-i? sources are all 
found to be steep, ax ~ 1-0). 

• In a few well-resolved cases (e.g. 3C 351, H02) the local value of R would be even higher in 
places than the integrated value we quote, representing an even greater challenge for SSC. 

• The fact that a synchrotron model can be fitted through the radio, optical and X-ray points 
in some high-i? sources would have to be a coincidence in an SSC model, although this is a 
weak constraint; the optical emission might also be SSC, as it is thought to be in a couple of 
low-R sources (Hardcastle 2001; Brunetti 2002). 

• An SSC model cannot explain the effects of beaming, if these are real; beaming suppresses 
SSC emission, so that we should see a weak anti-correlation with proxies of beaming like 
core prominence, at least for hotspots on the jet side. 

• The special status of i? = 1 suggests that there are few hotspots with B > Beq'. it is not 
obvious why the departures from equipartition should all be in the sense B < Beq. 

• There is no obvious mechanism that fully explains the observed luminosity dependence of 



A synchrotron model is a good fit to the overall spectrum of some of the most extreme sources, 
such as 3C 390.3 and 3C403; if we accept that X-ray synchrotron emission is possible in some 
hotspots, as it certainly is in the jets of FRl sources (e.g., Hardcastle et al. 2001b), then it may 
contribute to many of them. If the ability of a hotspot to produce X-ray synchrotron emission 
depended on its luminosity, then it might be the case that all hotspots have inverse-Compton emis- 
sion at a level consistent with R — 1 and equipartition magnetic fields, but that the low-luminosity 
hotspots have an additional synchrotron component that may greatly exceed the inverse-Compton 
emission. 



R. 



4.2. Luminosity-dependent synchrotron emission 
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It has already been argued in studies of optical synchrotron hotspots (e.g., Meisenheimer et 
al. 1997; Brunetti et al. 2003) that the high-frequency break in the synchrotron spectrum is a func- 
tion of hotspot luminosity, in the sense that optical emission is much commoner from hotspots 
of low radio luminosities. This fact can be explained (Brunetti et al. 2003) in terms of the lower 
synchrotron loss rates in the lower magnetic fields (assuming equipartition) and lower photon den- 
sities in hotspots of lower radio luminosity; we have already seen (Fig. 6) that there is a correlation 
between these quantities and R in our objects. In a standard hotspot spatial/spectral model (e.g.. 
Heavens & Meisenheimer 1987) the break in the synchrotron spectrum comes about when we aver- 
age over both the acceleration region itself and the regions downstream of it, in which synchrotron 
and inverse-Compton losses have had time to have an effect; we would not expect to see a break 
if we could resolve the acceleration region from the downstream emission. The high-frequency 
cutojf'm the synchrotron spectrum is a direct indicator of physics in the acceleration region, and 
results from inefficiency in particle acceleration at high energies: most importantly, from our point 
of view, particle acceleration will become inefficient if the energy loss timescale in the accelera- 
tion region (due to synchrotron and inverse-Compton losses) becomes shorter than the acceleration 
timescale. The fact that essentially all our X-ray data points fall on or below the line of an extrap- 
olation from the radio with ckrx = 10 (§3) shows that any luminosity-dependence of synchrotron 
radiation in the X-ray hotspots cannot simply be an effect of a changing frequency of the spectral 
break, as Brunetti et al. (2003) argued for optical hotspots; a = 1.0 is the canonical spectral index 
above the break, so that the break alone cannot produce the effect we see. This interpretation is 
supported by observations of well-studied luminous hotspots such as 3C405's, where the over- 
all spectrum requires the cutojf to be below the X-ray region. We conclude that the luminosity 
dependence of R must be an effect of the synchrotron cutoff. 

A full calculation of the frequency of the cutoff Uc depends on poorly known quantities such 
as the magnetic field strength in the acceleration region and the diffusion coefficient of relativistic 
particles. In the simplest case, with uniform magnetic field throughout the hotspot, a diffusion 
coefficient independent of both magnetic field and electron energy, and a non-relativistic shock, it 
can be shown that Vc B/(^B'^ + B^qY, with the constant of proportionality depending on the 
numerical value of the diffusion coefficient, where B is the magnetic field strength and Bic the 
equivalent inverse-Compton field strength, defined as Bic — ^/QJiqUic, with Ujc being the energy 
density in all photon fields. Pc in this calculation does exhibit a change over the parameters of 
the hotspots we have studied that would be sufficient in magnitude to explain the observed effect, 
although the detailed correlation with R is not particularly good (Fig. 9). Brunetti et al. (2003) 
consider two cases with less simplistic forms of the diffusion coefficient (Kolmogorov and Bohm 
diffusion coefficients), and are able to calculate corresponding numerical values for 7inax (their 
eq. 5). The functional form of 7niax for the Bohm coefficient means that the cutoff frequency is 
essentially constant, and (for the values quoted by Brunetti et al.) lies well above the X-ray region 
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for all our hotspots; the Kolmogorov coefficient gives a Uc lying well below X-ray frequencies for 
equipartition magnetic field strengths in all the hotspots. This illustrates the strong dependence of 
the expectation on the unknown microphysics of the acceleration process. The basic principle of 
this model remains plausible: hotspot luminosity (and therefore magnetic field and photon energy 
density) are controlling the high-energy cutoff of the synchrotron spectrum. 

Any relationship between R and core prominence is not easy to explain in a model where 
much of the X-ray emission is synchrotron. R does not have a simple dependence on beaming 
parameters in this model: the expected synchrotron flux of the hotspot increases with beaming, but 
so will the predicted inverse-Compton flux density (from both SSC and CMB scattering), since the 
prediction we make is based on the observed radio flux density and takes no account of beaming. 
We calculated the expected variation of R with angle to the line of sight 9 for a source whose 
intrinsic (rest-frame) properties were held constant. For modest beaming factors, corresponding 
to i'/c ~ 0.3, we find that it! does indeed increase as 9 gets smaller, but only for extremely low- 
luminosity hotspots, in which scattering of CMB photons is the dominant IC process; for more 
luminous hotspots the trend is reversed, and we would expect R to be largest for hotspots that are 
beamed away from us (that is, on the counterjet side of beamed sources). In any case, the amount 
of variation introduced by this process into the R value distribution is small, no more than a factor 
2 between minimum and maximum values, for v/c^ 0.3. Much higher speeds (v/c > 0.9) would 
be required, for reasonable hotspot luminosities, to obtain the order of magnitude scatter in the R 
parameter (after accounting for the luminosity dependence) that appears to be present in Fig. 5. 



4.3. Deceleration and beaming 

We have already shown (§3) that the apparent relationship between proxies of beaming, such 
as core prominence, and the X-ray brightness of the hotspot (H02; Georganopoulos & Kazanas 
2003) may at least partly be a selection effect in the available X-ray data. The current sample of 
X-ray hotspots, at least at low hotspot luminosities, is strongly biased towards beamed objects, 
while beamed objects are known to have brighter, more compact, flatter-spectrum hotspots on the 
jet side. It is now known that there are narrow-line radio sources that should lie close to the plane 
of the sky that have high R, and a few examples of sources (for example, 3C228 and 3C321) 
where hotspots on both sides of the source have high R. It is not clear, therefore, whether there is 
any beaming effect that needs to be explained from an X-ray perspective. 

If there is any beaming effect, then the two models discussed so far both have difficulty in 

explaining it, so alternative models must be considered. The standard way of explaining the (effec- 
tively) high R values in the X-ray jets of core-dominated quasars is to invoke highly relativistic bulk 
speeds and the consequent boost of the energy density of the CMB in the rest frame of the jet. We 
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regard this model as untenable in the present case, for several reasons. Firstly, large bulk Lorentz 
factors are required for even moderate R values (see the discussion of 3C 351 in H02) and these 
in turn constrain the source to lie at a small angle to the line of sight; this cannot possibly be the 
case for all or even most of our high-i? objects, which are drawn from a low-frequency-selected, 
lobe-dominated sample, and which include objects that, in unified models, must be close to the 
plane of the sky. Secondly, in the standard picture, the radio emission from the hotspots comes 
from the post-shock region, and so high bulk Lorentz factors are hard to achieve; although there 
are some effects that are best explained by moderate relativistic beaming in the post-shock flow, 
bulk Lorentz factors ~ 10 have never been required by observation, and are in fact inconsistent 
with the known properties of hotspots. 

A more viable model involving beaming effects is that of Georganopoulos & Kazanas (2003) 
(see §1). The picture they describe almost certainly has to be true at some level, but a quantitative 
test is difficult, since it relies on knowledge of the velocity and electron density structure of the jet 

upsteam of the hotspot that is hard to obtain observationally. In addition, this model cannot account 
for all the features of our data, such as the hotspot luminosity dependence of R. However, our data 
clearly do not rule out a beaming effect at some level and, if it is present, the standard inverse- 
Compton or synchrotron models cannot account for it without involving large speeds. A full test 
of this type of beaming model must await an unbiased sample of hotspots in which orientation and 
luminosity effects can be clearly separated. 



As discussed in §1, we have tried to distinguish between hotspots, defined as structures where 
the well-coUimated flow of the jet terminates, and jet knots, where the assumption is that the 
jet continues more or less unaffected by whatever process produces the increase in synchrotron 
emissivity. The key physical differences between the two systems are (1) that there is probably not 
a strong shock in FRU jet knots, since there is little evidence that the jets decelerate there, and (2) 
that the particles in jet knots probably have a shorter dwell time in the region of interest, since the 
downstream flow speed is likely to be faster, which could lead to spectral differences even if the 
acceleration processes are similar. In practice the distinction between the two types of feature is 
difficult to draw observationally: there are several features that we have considered to be hotspots 
in our sample (e.g., 3C 390.3 N, 3C 403 F6, 3C 275. 1 N) that might weU be jet knots in which the 
continuing jet is poorly defined. Equally, it must be the case that there is continued coUimated flow 
out of primary hotspots in cases where there is optical or possible X-ray synchrotron emission in 
the secondary hotspot, requiring in situ particle acceleration - 3C 35 1 's hotspots J and K are a good 
example. We see no observational differences between these borderline jet knot/hotspot sources 



5. Hotspots and jet knots 
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and clearly defined terminal hotspots. 

Should we therefore try to apply the results of the present work to jet knots as well as to 
hotspots? As we have argued above (§4.3) the generally favored jet X-ray emission mechanism for 
core-dominated quasars cannot apply to more than a small subset of our sources, and particularly 
not the narrow-line objects, some of which exhibit either possible jet-related X-ray knots (e.g. 
3C403) or clear jet-related X-ray emission (e.g. 3C321 and 3C452 in Appendix A). The X-ray 
emission mechanism here seems likely to be synchrotron, as in the jets of low-luminosity FRl 
sources. A full analysis of the known FRll jet-related X-ray emission is beyond the scope of 
the present paper, but from our work on hotspots we can make the 'prediction' (borne out by the 
obserations that we are currently aware of) that jet-related X-ray synchrotron emission in FRlls 
will be seen mostly in low-luminosity jet knots, and therefore should be particularly easy to find in 
low-luminosity FRII sources. The hotspot behaviour is also qualitatively similar to what is seen in 
the jets of some powerful quasars, such as 3C 273 (Sambruna et al. 2001; Marshall et al. 2001), in 
which the X-ray-to-radio ratio of jet knots decreases as the knot radio flux density increases. 



6. Conclusions 

We have shown that the properties of the X-ray emission of hotspots depend strongly on their 
overall radio luminosity. High-luminosity hotspots, of the type originally examined in inverse- 
Compton studies, consistently show X-ray emission that is close to being consistent with the 
predictions of a synchrotron self-Compton model with an equipartition magnetic field. Low- 
luminosity hotspots sometimes (and maybe always) have X-ray emission that is much brighter 
than would be expected in this model. We argue that: 

• The good agreement between IC models and data seen for the luminous hotspots continues 
to suggest that the X-ray emission mechanism in these systems really is synchrotron self- 
Compton, that magnetic fields really are in equipartition, and that populations of protons 
which dominate energetically by large factors (> 100) and/or very small filling factors are 
not present. 

• Models in which the unexpectedly strong X-ray emission from some low-luminosity hotspots 
indicate a large departure from equipartition in these objects are not plausible for a number 
of reasons: synchrotron emission is more likely. 

• If a synchrotron model is adopted, the high-frequency cutoff of the synchrotron spectrum 
must be dependent on luminosity in order to explain the X-ray emission from all hotspots. 
This is physically plausible, but a fully quantitative test depends on the microphysics of the 
acceleration process. 
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• There is little significant evidence that relativistic beaming is important in the current sam- 
ple; an unbiased sample of X-ray hotspots would be of great importance in testing beaming 
models. 

• It may be possible to extend our conclusions on hotspots to the jet-related X-ray features 
seen in a number of FRII sources, particularly those at relatively large angles to the line 
of sight; if so, we would expect that they would show the same luminosity dependence, in 
the sense that only low-luminosity jets would show strong X-ray synchrotron emission. A 
synchrotron origin for the jets in these sources would suggest a continuity between their 
properties and those of the lower-power FRIs, for which a synchrotron interpretation is well 
established (e.g., Hardcastle et al. 2001b). 
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Fig. 1. — Distribution of the ratio R of observed X-ray flux to IC prediction for the Chandra 
sample. Left: distribution for X-ray detected hotspots. Right: distribution for compact hotspots 
with no X-ray detection (sources can move to the left). 
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Fig. 2. — B! is plotted against R, where R is the ratio between the observed X-ray flux density and 
the prediction of an inverse-Compton model at equipartition, and R' is defined as the ratio of the 
observed flux to the extrapolation of the radio flux density assuming a = 1.0, i.e. to the amount 
of X-rays that could (conservatively) have been produced by synchrotron emission with a straight 
spectrum. Diagonal arrows show upper limits from non-detected X-ray sources. Insets show the 
broad-band SEDs of two extreme sources on the plot: top left, Cygnus A hotspot A; bottom right, 
3C 390.3, N hotspot. The data points are from the literature or from maps available to us, the solid 
line represents the best-fitting synchrotron model, and the dot-dashed and dotted lines represent 
the equipartition synchrotron self-Compton and CMB inverse-Compton models, respectively. 
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Fig. 3. — Spectral index difference (curvature indicator) against R for sources with optical hotspots 
or upper limits on optical flux. The error bars show the statistical errors on X-ray flux density only, 
as these are the dominant errors. 
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Fig. 4. — it! plotted against the total rest-frame 178-MHz luminosity, from the 3C/3CRR measure- 
ments. Since most sources have more than one hotspot, typically two R values are plotted for a 
given source luminosity. 
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Fig. 5. — it! plotted against the rest-frame 5-GHz luminosity of the hotspot. The dotted lines 
extending down from the data points show the approximate lowest value of R that could have been 
detected with the data, assuming a nominal Chandra sensitivity of 0.1 nJy at 1 keV. 
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Fig. 6. — R plotted against the magnetic field and photon energy densities in the hotspots. The 
plotted photon energy density takes into account both synchrotron and CMB photons. Note that 
the two plots have different scales on the x-axis. 
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Fig. 7. — R plotted against core prominence. Stars indicate broad-line objects (broad-line radio 
galaxies and quasars) that are expected to lie at angles < 45° to the line of sight in unified models. 
Circles around data points indicate hotspots on the same side of the source as a distinct one-sided 
radio jet. 
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Fig. 8. — Core prominence plotted against hotspot luminosity for the X-ray hotspot sample. 
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Fig. 9. — R plotted against the synchrotron cutoff frequency for the hotspots in the sample, us- 
ing the proportionality between Uc and magnetic field strength for a constant diffusion coefficient 
quoted in the text. The constant of proportionality (i.e. the normalization of the X-axis) is chosen 
simply to illustrate that the magnitude of the effect could be significant in this situation, and has no 
physical basis. The magnetic field strength used is the equipartition field derived from our models, 
and the photon energy density is a combination of the microwave background and the integrated 
synchrotron spectrum for each hotspot. A trend is apparent, in the sense that low-i? hotspots have 
low cutoffs while high-i? ones have high cutoffs, but the scatter is large. 
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Table 1. 3C FRII sources observed with Chandra 
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Table 1 — Continued 



Source 
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Note. — 5'i78 is the 178-MHz flux density, on the scale of Baars et al. (1977), mostly taken from 
Laing et al. (1983) or Spinrad et al. (1985). ckr is the low-frequency spectral index, typically between 
178 and 750 MHz. Types are based on optical and emission-line characteristics, and are as follows: E, 
low-excitation radio galaxy; N, narrow-line radio galaxy; B, broad-line radio galaxy; Q, quasar. Values 
of S'corc,55 the 5-GHz core flux density, are mostly taken from the compilation on the 3CRR web pages 
(http://www.3crr.dyndns.org/) or from the radio maps referred to in this paper. All Chandra datasets are 
from the ACIS-S except 3C295, where the ACIS-I was used. No grating data met our selection criteria. 
Livetimes are the filtered times if filtering was carried out, and the uncorrected livetime otherwise. 
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Table 2. VLA radio observations used in this paper 



Source 


Proposal ID 


Frequency 


Time on 


Date 


Reference 






(GHz) 


source (h) 




(if published) 


i.c a 1 


AVXQn 
AtdoKJ 


o.D 


i 


1 QQQ Aiirr no 

lyyy Aug uz 








o.D 


lb 
i 


'~)(\(\(\ Ton 1 fi 

zuuu J an i o 




y 


A 13 "2 AO 

Aojoy 


A Q 

4.y 


A < 

4.D 


1 QfiA f\/l 

lyoo May U4 


D 


An 


Ad /yo 


/I fi 


1 


lyyo iNov u/ 




1 no 




0.4 






i 






Q /I 
0.4 






z 


o/^ 1 1 




Q /I 
0.4 






z 


1 /y 




/I Q 

4.y 


U.D 


lyoD iviar zi 








A Q 

4.y 






-2 
D 


ZUU 




fi A 
0.4 






1 
1 




Ad /yo 


fi < 
O.D 


-3 
D 


1 OOA M/-v^r fifi 

lyyo iNov uo 






Ad /yo 


fi < 
O.D 


A 

4 


lyyo iNov u/ 








A fi 
4.0 






4,D 


jL- ziy 




L.J 






A 
4 


jL- ZZU. i 




fi A 
0.4 






O 


Ji^ ZZo 




fi /I 
0.4 






1 
i 


'XC T^/l 
jL- Zj4 




A Q 

4.y 


U.D 


1 QfiO tToU A1 

lysy rcD ui 




ZOj 




4.y 






D 


JK^ ZOD 


Ar ioO 


A fi 
4.0 


-3 
D 


1 QOn A t->i- 00 

lyyu Apr zz 


1 1 
1 1 


'ir' otn 1 




fi A 
0.4 






1 
1 


3C280 


AVI 57 


8.4 


0.7 


1988 Dec 22 




3C281 


AB631 


1.4 


0.3 


1992 Nov 18 




3C294 


AM224 


4.7 


3.5 


1987 Oct 11 




3C295 




8.4 






1 


3C303 


KRON 


4.9 


0.3^ 


1981 Apr 20 


10 


3C321 


AV127 


4.8 


3.7 


1986 Apr 10 








1.5, 4.8 


2.9, 3.7^ 


1986 Aug 29 




3C324 


AF186 


4.9 


3 


1990 Apr 22 


11 


3C330 




8.4 






1 


3C334 




4.9 






5 
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Table 2 — Continued 



oouice 


-Tioposai HJ rTcCjuency 
(GHz) 


iinie on uaie Kcierence 
source (h) (if published) 


3C351 


8.4 


1 


3C 390.3 


1.5 


4 


3C401 


8.4 


2 


3C403 


8.4 


7 


3C405 


4.5 


8 


3C 427.1 


8.4 


1 


3C438 


8.4 


2 


3C452 


8.4 


7 


Pictor A 


4.9 


9 



Note. — We list VLA observational details only for observations that we have retrieved 
from the archive and reduced ourselves in the course of this project; for other observations 

we were able to obtain electronic maps from others (or already had them ourselves) and 
the reader is referred to the references given below for the observational information. All 
data retrieved from the archive were taken with the VLA in its A configuration, except 
where otherwise noted. 

References. — (1) Gilbert et al. 2004; (2) Hardcastle et al. 1997; (3) Belsole et al. 2004; 
(4) Leahy et al. 1998 (the 3CRR Atlas); (5) Bridle et al. 1994; (6) Worrall et al. 2001; (7) 
Black et al. 1992; (8) Perley, Dreher & Cowan 1984; (9) Perley, Roser & Meisenheimer 
1997; (10) Kronberg 1986; (11) Femini et al. 1993. 

*Only one observing frequency of 12.5 MHz bandwidth was used. 

''B -configuration data. 
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Table 3. Radio and X-ray flux densities and predicted inverse-Compton flux densities 





l-rr\f CT\r\f 


A n (T CI 


j-vjnz, iaciivj 


1 -IfpV flllY 


Pff^Hi /^f f^H fliiY 

-TiCUlV^LCLl ilUA 


Ratio /? 
ivaLiVj XL 








flllY 

iiLiA \jy ) 






yyjUcid. VCU./ 

trr^^ri i r*ff^r\ \ 
jji CU-IC LCU- J 


3C6.1 


N 


0.36 


0.340 


0.45 


0.19 


2.3 




s 


0.41 


0.200 


0.09 


0.081 


1.1 


3C9 


N 


0.38 


0.038 


< 0.09 


0.024 


< 3.6 




s 


0.39 


0.012 


< 0.09 


0.0046 


< 20 


3C47 


s 


0.434 


0.181 


0.54 


0.040 


14 




N 


1.89 


0.127 


< 0.1 


0.015 


< 6.6 


3C109 


s 


0.377 


0.181 


0.15 


0.033 


4.6 




N 


0.274 


0.007 


< 0.09 


0.00053 


< 169 


3C123 


E 


1.1 X 0.54 


5.12 


4.6 


2.2 


2.1 




W 


1.0 X 0.13 


0.341 


0.18 


0.059 


3.1 


3C 173.1 


s 


0.83 


0.033 


0.2 


0.0022 


91 




N 


0.26 


0.009 


< 0.12 


0.00051 


< 237 


3C179 


W 


0.145 


0.063 


1.54 


0.026 


63 




E 


0.45 


0.038 


< 0.26 


0.0070 


< 37 


3C200 


N 


0.6 


0.057 


< 0.1 


0.011 


< 9.5 


3C207 


E 


0.27 


0.044 


0.69 


0.0095 


73 


3C212 


N 


0.144 


0.035 


< 0.14 


0.014 


< 10 




S 


0.25 


0.110 


< 0.14 


0.056 


< 2.5 


3C215 


E 


1.0 


0.012 


< 0.04 


0.0011 


< 36 


3C 220.1 


E 


0.27 


0.021 


< 0.13 


0.0025 


< 52 




W 


0.27 


0.023 


< 0.13 


0.0030 


< 44 


3C228 


N 


0.203 


0.070 


0.45 


0.019 


24 




S 


0.265 


0.132 


1.3 


0.042 


31 


3C254 


W 


0.29 


0.146 


0.54 


0.061 


8.8 


3C263 


E 


0.39 


0.582 


1.0 


0.25 


4.0 




W 


0.18 


0.023 


< 0.06 


0.0054 


< 11 


3C265 


E 


0.356 


0.272 


0.35 


0.16 


2.2 




W 


0.73 


0.048 


0.13 


0.088 


15 


3C 275.1 


N 


1.4 X 0.2 


0.191 


1.78 


0.093 


19 
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Table 3 — Continued 



Soiircp 


Mot snot 


AricT sizp 
( arcspc^ 


S-fiFf/ radio 
flux Qv") 


1 -keV flux 


Prpdictpd flux 


Ratio B 
DTpdictpd^ 




S 


0.378 


0.111 


< 0.12 


0.038 


< 3.1 


3C280 


E 


0.186 


0.082 


0.31 


0.046 


6.7 




W 


0.146 


0.631 


0.6 


1.2 


0.48 




Wc 


0.081 


0.035 


0.07 


0.021 


3.35 


3C281 


N 


1.04 


0.129 


< 0.16 


0.022 


< 7.3 


3C294 


N 


0.283 


0.143 


0.12 


0.17 


0.72 




S 


0.43 


0.022 


0.14 


0.0088 


16 


3C295 


N 


0.1 


1.29 


1.4 


0.78 


1.8 




S 


0.1 


0.92 


0.94 


0.85 


1.1 


3C303 


w 


1.1 X 0.28 


0.257 


4.0 


0.026 


154 




E 


0.57 


0.0025 


< 0.16 


0.00039 


< 408 


3C321 


E 


0.69 


0.125 


0.3 


0.006 


48 




W 


2.7 X 0.45 


0.020 


0.12 


0.00057 


210 


3C324 


E 


0.365 


0.277 


0.20 


0.21 


0.93 




W 


0.301 


0.085 


0.16 


0.040 


4.04 


3C330 


N 


0.45 


0.625 


0.35 


0.42 


0.81 




S 


0.20 


0.102 


0.068 


0.028 


2.4 


3C334 


S 


1.34 X 0.3 


0.018 


0.54 


0.0018 


292 




N 


0.5 


0.007 


< 0.4 


0.00066 


< 604 


3C351 


J 


0.16 


0.167 


4.3 


0.051 


85 




K 


0.8 


0.406 


3.4 


0.087 


39 




S 


0.16 


0.0025 


< 0.05 


0.00014 


< 362 


3C 390.3 


N 


1.3 X 0.5 


0.087 


4.5 


0.003 


1380 




S 


3.7 


1.30 


1.3 


0.07 


19 


3C403 


Fl 


0.275 


0.021 


1.0 


0.00047 


2149 




F6 


0.256 


0.035 


1.8 


0.0013 


1414 


3C405 


A 


2.5 X 1.25 


38.0 


19.4 


15.5 


1.3 




B 


0.44 


3.04 


4.5 


0.71 


6.3 




D 


1.09 


30.3 


29.2 


13.6 


2.1 
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Table 3 — Continued 



Source 


Hotspot 


Ang. size 


5-GHz radio 


1-keV flux 


Predicted flux 


Ratio R 






(arcsec) 


flux (Jy) 


(nJy) 


(nJy) 


(observed/ 














predicted) 




E 


0.45 X 0.63 


1.68 


1.2 


0.23 


5.2 


3C 427.1 


N 


0.165 


0.019 


< 0.17 


0.0030 


< 56 




S 


0.14 


0.025 


< 0.17 


0.0046 


< 37 


3C452 


W 


0.705 


0.033 


0.34 


0.00095 


356 




E 


3.0 


0.067 


< 0.05 


0.0039 


< 13 


Pic A 


W 


0.75 


1.93 


89 


0.20 


454 




E 


0.75 


0.467 


< 0.16 


0.044 


< 6.4 



Note. — The hotspot identifier is usually N, S, E or W, referring to the obvious or brightest 
hotspot in the north, south, east or west lobes. Exceptions are made where a multiple-hotspot 
source has names for the individual components that are used relatively widely in the literature; 
this is true of 3C 405 (notation of Hargrave & Ryle 1974), 3C 351 (notation of Bridle et al. 1994) 
and 3C403 (notation of Black et al. 1992). The angular sizes quoted are the radii of homoge- 
neous sphere models fitted to the radio data, as described in the text, except where two numbers 
are quoted, in which case they are the length and radius of a cylinder and are generally directly 
measured from high-resolution maps. The measured hotspot radio flux densities have been scaled 
to a lab-frame radio frequency of 5 GHz using a spectral index a — 0.5 for ease of comparison. 
The 1-keV flux densities are the values inferred from spectral fitting or the observed count rate, as 
described in the text, and are the unabsorbed fluxes (assuming Galactic absorption). 
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Table 4. Optical flux densities used in this paper and associated HST observational details 



Source 


HS 


Freq. 
Hz) 


Flux 
density 

(/^Jy) 


Origin 


Ref. 


Obsid 


Filter 


Time on 

source 
(s) 


Date 
observed 


'ic An 


e 


J.D 


< U.o i.UJ 


JJCT 


1 


U44yZlUl 


rjjj W 


OUU 


1 OOO Ton 

lyyy jan ju 


ji_ iuy 


e 


A K 




JJCT 


1 


Tn7T 1 Qr»i 


r /UZW 


JOU 


1 QQ< Aiirr lA 

lyyj Aug Z4 




11 


c c 
J.J 


< Z.O 


JJCT 


1 
i 


U44y4oUl 


rjjj W 


OUU 


1 QQQ A r^r n'C 

lyyy Apr uj 




W 




< U.40 (^-C i.Oj 














'?r' 1 1^ 1 
j^-- i / J. i 


c 
o 


A 'X 


< i.Z -'-•4J 


tlol 


/ 


TT07T OOni 
UZ /LZ^JUl 


r / UZ W 


JUU 


1 QQ/I Till 07 

iyy4 jui z / 




W 


c c 
J.J 


< u.40 ^< <J.o 1 ) 




/ 


U44yjL-Ul 


r J J J w 


DUU 


1 QQQ A/Tqi- n^^ 
lyyy ividr uo 




C 


J.J 


U.OO l^<. U.4Dj 


nol 


7 

/ 


U44yo /Ul 


r J J J VV 


DUU 


1 QQQ Ton 1 S 
Lyyy Jdn 1 o 




IN 


'X ^ 
J.J 


^ n 8P; n Q/l^ 


riijl 


/ 


U Dr'A J / U i 


1h78^T P 
r' / oji^r 


zuuu 


onm \/fQA7 
zuui iviay ju 




e 




i.U4 (^i.i4J 
















w 

W 


c c 
J.J 


< u.oo U.4iJ 


iliji 


7 


U44yUvyUi 


r'jjj W 


Ann 

OUU 


1 QQQ A^^of 1 fi 

lyyy iviar lo 




a 


A K 
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U.O (,U.Oj 


TJCT 


c 
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UZoilUZUi 


ro/j W 


lUUU 


1 QQ/^ ThoK 1 8 

lyyo reo lo 
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«: « 
J.J 
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UZV^ 1 UJUZ 
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^ 7nn 
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T QQ*;; Anr 01 
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< i.U ^< -L- -LJ 




7 


TT07T /Ipm 
UZ /L4rUl 


r / UZ W 


JUU 


1 QQ^ A/fa\7 1 

lyyj iviay iz 


Q/^ 07^ 1 
Jl-- Z / J. 1 


IN 


A 7 


U.44 (U.4o ) 




7 


UZoiiUjUl 


T7^i7 ^\X/ 

ro/JW 


ioUU 


1 QQ^ Till 0*; 

lyyj Jul zj 


ZoU 




A Q 


n 'JO ((\ 'XA\ 
U.JZ (U.j4; 




7 


UZvjAUoUi 


rOZZ w 


ooUU 


1 QQ/l A iirr 00 

iyy4 Aug zz 




W 
w 




U.4 (^<^ U.4J 
















wc 




< U.4 U.4J 














'XC OQ'n 

zyj 


M 
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A 'X 
4. J 


U.U/O yj.KJOL) 




/I 
4 


UZC4UAU1 


r /UZW 


iZOUU 


1 QQA Ton 1 A 

lyyo Jan 14 




c 




U.UZ \ \j.\jL) 














3C303 


w 


5.5 


7.5 


OHP 


2 










3C330 


N 


5.5 


< 0.5 (< 0.6) 


HST 


5 


U3A14X01 


F555W 


600 


1996 Jun 03 




S 




< 0.5 (< 0.6) 














3C334 


s 


5.5 


< 0.7 (< 0.86) 


HST 


7 


U4492V01 


F555W 


600 


1998 Dec 27 


3C351 


J 


4.3 


2.4 (2.6) 


HST 


5 


U2X30601 


F702W 


2400 


1995 Nov 30 




K 




1.9 (2.1) 














3C 390.3 


N 


4.5 


2.2 (2.6) 


NOT 


6 








1996 Jul 18 


3C403 


Fl 


4.3 


0.66(1.08) 


HST 


8 


U27L7601 


F702W 


280 


1994 Jun 26 




F6 




1.32 (2.16) 














3C405 


A 


4.6 


< 80 


Calar 


2 








1985 Oct 
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Table 4 — Continued 



Source 


HS 


Freq. 

(xlO^^ 
Hz) 


Flux 
density 

(/^Jy) 


Origin 


Ref. 


Obsid Filter Time on 
source 
(s) 


Date 
observed 




B 




< 46 


Alto 


2 








D 




< 5 




2 






Pictor A 


W 


4.5 


130 


ESO 3.6m 


3 




1985 Nov 08 



Note. — Hotspots are identified as in Table 3. Optical flux densities are corrected for Galactic extinction 
if only one value is given; where two are given the second (in parentheses) is the corrected value. HST 

observational details are given where HST data were used; other data points are taken from ground-based 
observations described in the literature. Data points for the same object have the same origin unless dif- 
ferent origins are explicitly listed in the Table. NOT indicates the Nordic Optical Telescope and OHP the 
Observatoire de Haute Provence. 

References. — (1) Hardcastle et al. 2001a; (2) Meisenheimer et al. 1997; (3) Roser & Meisenheimer 
1987; (4) Harris et al. 2000; (5) H02 (6) Harris et al. 1998; (7) This paper; (8) R. Kraft et al. , in prep. 
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A. Newly detected hotspots 

Below we present images of the hotspots newly detected in the course of this work that are 
not expected to be discussed in more detail in other papers. 

A.I. 3C6.1 

The N hotspot of this narrow-line source is the clearest detection (Fig. 10), but there is a weak 
detection of the S hotspot too, at well over 3a significance. Some X-ray emission is associated 
with the lobes. 



A.2. 3C47 

This quasar's bright S hotspot (on the jet side) is detected (Fig. 11), but there is no obvious 
detection of the fainter N hotspot. Extended emission is clearly visible in the X-ray image: since 
it is extended in the direction of the lobes some of it may well be inverse-Compton emission, but 
probably a large fraction of it comes from a cluster environment, particularly as 3C 47 exhibits a 
strong Laing-Garrington effect (e.g., Leahy 1996). 



A.3. 3C 109 

The southern hotspot of this broad-line radio galaxy is detected; the southern side is the jet 
side and a weak jet can be traced into the hotspot (Gilbert et al. 2004). Some excess extended 
emission from the lobes can be seen in these images, and is consistent with inverse-Compton 
emission at approximately the level expected from equipartition in the lobes. 

A.4. 3C 173.1 

There is a weak detection of the S hotspot of this low-excitation radio galaxy (Fig. 13), at 
well over 3a significance. The extended emission here is again probably a combination of inverse- 
Compton emission and a thermal environment. The hotspot here is on the counterjet side. 
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Fig. 10. — The X-ray hotspots of 3C6.1. The greyscale shows the 0.5-5 keV Chandra counts 
smoothed with an O'.'S FWHM Gaussian; black corresponds to 1 count per 0'.'246 pixel. The 
contours are of the 8.4-GHz B-configuration VLA map at 0'.'94 x 0'.'53 resolution, and are at 
0.2 X (1,4, 16. . .) mJybeam-i. 
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Fig. 11. — The southern X-ray hotspot of 3C47. The main greyscale (left) shows the 0.5-5 keV 
Chandra counts smoothed with an Qf!5 FWHM Gaussian; black corresponds to 5 counts. The 
contours are of a I'.'O resolution 1.6-GHz VLA map taken from Leahy et al. (1998), and are at 
0.3 X (1, 4, 16 . . .) mJy beam~^. The inset (right) is the same X-ray image with the same greyscale 
level, but with contours from the 0'/39 x 0'/36 resolution 4.8-GHz VLA map at 80 x (1, 4, 16 . . .) 
liJy beam~^. 
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Fig. 12. — The southern X-ray hotspot of 3C 109. The main greyscale (right) shows the 0.5-5 
keV Chandra counts smoothed with an O'.'S FWHM Gaussian; black corresponds to 1.5 counts per 
0'.'492 pixel. The line across the image is the Chandra readout streak. The contours are of a 2'.'5 
resolution 8.4-GHz VLA map taken from Gilbert et al. (2004), and are at 0.2 x (1, 4, 16 . . .) mJy 
beam~^. The inset on the left shows the hotspot; the X-ray map is the same but contours are from 
a 0'.'25 resolution 8.4-GHz map also from Gilbert et al. (2004), at 0.1 x (1, 4, 16 . . .) mJy beam~\ 
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Fig. 13. — The southern X-ray hotspot of 3C 173.1. The greyscale shows the 0.5-5 keV Chandra 
counts smoothed with an I'.'O FWHM Gaussian; black corresponds to 4 counts per 0'.'246 pixel. 
The contours are of a 1'.'7 resolution 8.4-GHz VLA map taken from Hardcastle et al. (1997), and 
are at 0.2 x (1, 4, 16 . . .) mJy beam"\ 
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A.5. 3C321 

Both hotspots of this nearby narrow-line FRII source are detected, weakly but convincingly 
(Fig. 14). In addition, there is X-ray emission from a weak radio jet entering the S hotspot from 
slightly W of N (just visible on our high-resolution contour map). The nuclear region shows very 
unusual structure. The component associated with the radio core is extended, and, if the relative 
astrometry of the Chandra and radio data is correct, shows quite strong X-ray emission from 
the radio-weak jet pointing SE; the NW compact bright X-ray source is positionally coincident 
not with the bright NW radio jet, which is not clearly detected in X-rays, but with the nearby 
companion galaxy seen with HST (Martel et al. 1999) if we align the center of the host galaxy 
with the radio core and the brightest component in the X-ray (the HST data has the usual arcsec- 
scale astrometric offset, and there are no obvious independent features with which to align the 
two datasets). There is also extended emission around the two galaxies that appears to be spatially 
coincident with the known optical line-emitting material (Baum et al. 1988). The strongly different 
jet-counterjet asymmetry in the radio and X-ray is hard to explain in a model in which the jet and 
counterjet are intrinsically symmetrical. 

A.6. 3C324 

Both hotspots of this small narrow-line radio galaxy are detected (Fig. 15), the E hotspot 
clearly, the W one more marginally. 

A.7. 3C452 

The W hotspot of this low-redshift narrow-line radio galaxy is clearly detected (Fig. 16), as 
is a faint linear X-ray feature pointing W from the nucleus towards the known radio jet in the W 
lobe (though the features shown in Fig. 16 have no detected radio counterparts on high-resolution 
maps). Extended emission associated with the lobes has already been reported and is modeled in 
terms of inverse-Compton emission (Isobe et al. 2002). 
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Fig. 14. — The X-ray hotspots of 3C321. The main greyscale shows the 0.5-5 keV Chandra 
counts smoothed with an 2'.'0 FWHM Gaussian; black corresponds to 1 count per 0'.'492 pixel. The 
contours are of a 1.4-GHz VLA map with 15" x 13" resolution, and are at 2 x (1, 4, 16 . . .) mJy 
beam~^ Insets show the same map smoothed with an 0"5 FWHM Gaussian, and contours from a 
4.8-GHz VLA map with 0"45 x 0"40 resolution, at 0.15 x (1, 4, 16 . . .) mJy beam"^; black is 1 
count per pixel for the hotspots and 5 counts per pixel for the nuclear inset. 
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Fig. 15. — The X-ray hotspots of 3C324. The greyscale shows the 0.5-5 keV Chandra counts 
smoothed with an Qf!5 FWHM Gaussian; black corresponds to 4 counts per (y/246 pixel. The 
contours are of the 4.8-GHz VLA A-configuration map at (X'SQ resolution, and are at 0.15 x 
(1,4, 16...)mJybeam-^ 
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Fig. 16. — The western X-ray hotspot of 3C 452. The main greyscale shows the 0.5-5 keV Chandra 
counts smoothed with an 2'.'0 FWHM Gaussian; black corresponds to 10 counts per 0'.'246 pixel. 
The contours are of an 8.4-GHz VLA map with 2'.'5 resolution from Black et al. (1992), and are at 
0.3x(l,4,16...) mJy beam~^. Insets show the same X-ray image smoothed with a O'.'S FWHM 
Gaussian, with black being 1 count per 0^246 pixel. Right, the hotspot: contours are of an 8.4-GHz 
VLA map with 0'.'25 resolution from Black et al. (1992), and are at 0.1 x (1, 4, 16 . . .) mJy beam~^ 
Top, the inner jet (no contours are shown). 
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B. Synchrotron and inverse-Compton emission 

To guide the reader in interpreting the physics of synchrotron and inverse-Compton emission, 
we include here a brief sketch of the underlying physics. In practice these calculations are carried 
out by the computer code discussed in the paper and in Hardcastle et al. (1998), but it is useful to 
set out the theoretical underpinning of the code's results and to outline the key dependences of the 
model parameters. 

We take as a fiducial assumption (which can then be tested by observation) the equipartition 
of energy between electrons (or more generally particles of all kinds) and magnetic field. If the 
electron energy spectrum (number per unit energy per unit volume) is described by a function 
N{E), then equipartition implies (in SI units) 

— = / EN{E)dE + u^^ (Bl) 

where Emin and i^max give the range of electron energies, B is the magnetic field strength, /xq 
is the permeability of free space and mnr is the energy density in non-radiating particles. It is 
conventional to let k be the ratio of the energy densities in non-radiating and radiating particles: 
then 

= (1 + / EN{E)dE (B2) 

It is easy to see that the (1 + k) term can also be used to describe an arbitrary departure from 
equipartition between the electrons and magnetic field. Our fiducial assumption is equivalent to 

K = 0. 

Now let us consider for simplicity a power-law distribution of electron energies, N{E)dE — 
NoE'PdE. Then the integral can be carried out analytically: 

S2 



where 



^{1 + k)NoI (B3) 



2-p 



ln(£;niax/^min) P = 2 

r;(2-p) _ ^(2-p)l 



In practice, as described in the text, we may use more complicated electron energy spectra, and 
then it is easiest to determine / numerically. 

The volume synchrotron emissivity of the ensemble of electrons at a given source-frame fre- 
quency u may be written (e.g., Longair 1994) 



(B4) 
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Here nie is the mass of the electron, e is its charge and c is the speed of Ught; eo is the permittivity 
of free space. 9 is the pitch angle of the electrons with respect to the magnetic field direction and 
X is defined by 



X — 



3e E^B sin 9 

F{x) is a sharply peaked function of x, reflecting the fact that electrons of a given energy radiate 
at a well-defined frequency; 

poo 

F{x) ^x Kr,/3iz)dz (B5) 

J X 

where 7^5/3 is the modified Bessel function of order 5/3. Assuming pitch angle isotropy, we can 
integrate equation B4 over pitch angle and our assumed electron power law (e.g., Longair 1994) to 
find that 

(p-i) (p+i) 

J{u) = CNqu — —B— (B6) 



where 

C = c(p)- 

eocme V e 

and c{p) is of order 0.05 and depends only weakly on p. Since equation B6 describes a power 
law in frequency, the electron energy index p can be determined by observation: typical values lie 
in the range 2-3, and our assumption of a low-frequency spectral index a — 0.5 corresponds to 

If we know the emissivity J, we can use equation B6 to eliminate A'^o from equation B2 

l^^{^ + i^)^^~B^I (B7) 

and we can now solve for B: 



B 



2 

P+5 

(B8) 



In practice, we take account of the fact that J(z/) is not always a power law by numerically inte- 
grating equation B4 and then solving equation B7 numerically with a root-finding algorithm, but 
the main dependences are encapsulated in equation B8. Since the energy density is proportional to 
B'^, we can see that it increases as (1 + thus, a non-zero value of k affects the magnetic 

field strength in the expected sense. Moreover, we can now substitute back into equation B6 to 
eliminate B: this gives 

J(i/)i^ = C^Nqu'^^^) [2//o(1 + k)/]^ (B9) 
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and, since J(i/) and v are known and constant for a given observation and a known source ge- 
ometry, we can see that the number density of electrons is expected to decrease with increasing 
k: 

Nq(x{1 + k)~^ (BIO) 

If the geometry is doubtful, the calculated emissivity is a function of volume: J oc S/V, where S 
is the observed flux density. So we expect 

No(xV~^ (Bll) 



For sphere of radius r with uniform particle and magnetic field density, the synchrotron self- 
Compton emissivity at a given frequency i/i is given by (Hardcastle et al. 1998) 

Jic(^i) = ^m.^cViaxr f^^^ ^^^^^^ f {x)dvodE (B12) 

where rrif. is the electron mass, c is the speed of light, ut is the Thomson cross-section, J{vq) 
is the synchrotron emissivity as a function of frequency, and f{x) is a function of E, ui and uq 
defined by Rybicki & Lightman (1979). The code we use performs this integration numerically 
for a given synchrotron and electron energy spectrum, and also carries out the similar calculation 
for illumination from the microwave background radiation; an analytical form of the integral for 
power-law electron and photon distributions could be derived but is not necessary here. The key 
feature of this equation is that the inverse-Compton emissivity (which determines the predicted 
IC flux density and thus R for a given observed flux) is linear in the number density of electrons 
N{E), and thus linear in A^o for the power-law analysis we have described above. For a given 
source, with known spatial and spectral properties, it is the dependence of A^o on (1 + k) given by 
equation BIO that primarily determines the value of R. (If the spectrum is not a pure power law, 
the change in the form of 7(1/0) as a result of the change in B also has a non-negligible effect.) 

If the volume is not known, for example because of a low filling factor, then the dependence 
of Jic on volume can be determined from above: 

Jic oc V^V'^V-^ = (B13) 

and this means that the observed inverse-Compton flux density, which is proportional to VJic, goes 
as T/^/3-V(p+5) This is a very weak dependence for plausible p values: for p — 2, Sic ^ 
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